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I  .  ^  Abstract 

;  Z'  f  1  .  ~ 

yc  derive  the  nost  general  e<|uatlons  of  notion  for  the  electrons  and  the  electronagnet Ic 
tii'ld  111“  tree  electron  laser  including  the  effects  of  diffractlor  and  pulse  propagation. 
The  t  i e ld''TrvTrtxitTon  I^expressed  In  terns  of  the  anplltudes  and  phases  of  a  complete  set  of 
transverse  nodes.  The  analytic  solution  is  given  In  the  small  signal  regime,  where  the 
thtMiry  Is  sliown  to  be  in  excellent  agreenent  with  a  recent  experinent , at  Ursay. 

Introduction  \ 

we  summarize  <  new  approach^  to  calculating  the  three  dimensional  effects 
electron  lasers.  The  previously  mentioned  approaches  consider  the  growth 
aTSng~tfte~ propagation  or  z  axis  by  evaluating  Its  change  at  each  point 
and  integrating  numerically  through  the  Interaction  region  In  the  time  domain^*'  or 
treijunncy  domain.  These  techniques  all  demand  long  conputer  runs  if  they  are  to  be 
applied  to  a  real  experimental  situation.  Our  approach  decomposes  the  problem  into  the 
**^  luninum  number  of  physically  observable  quantities:  the  transverse  optical  modes  of  the 
-.ystem.  The  field  evolution  is  expressed  in  terms  of  a  complete  set  of  orthogonal  trans¬ 
verse  nodes;  equations  are  developed  for  the  propagation  of  the  amplitude  and  phase  of  each 
node.  In  physical  systems  which  operate  on  a  few  of  the  lowest  order  modes,  this  approach 
;,reatly  increases  the  accuracy,  and  may  reduce  the  required  computer  time  for  the  calcula¬ 
tion  by  working  in  a  vector  space  well  matched  to  the  solution  ot  the  problem.  For  the 
oscillator  case,  the  appropriate  choice  of  nodes  Is  the  sot  of  eigennodes  of  the  cavity, 
lor  •  :n  1 1 1  ter ,  the  vector  space  of  nodes  is  deterinlned  by  the  characteristics  of  the 
ini.iit  nou»;,  which  IS  presumably  a  TEM  Gaussian  mode.  In  either  device,  an  optimum  design 
would  result  in  the  excitation  ot  as  V^w  of  the  higher  order  modes  as  possible.  The  nodal 
decomposition  method  is  therefore  well  adapted  to  the  prediction  and  optimization  of  the 
opiTution  of  the  tree  electron  laser. 

He  begin  hy  developing  the  most  general  equations  for  the 'propagation  of  the  mode  ampli¬ 
tudes  and  phases  in  a  FFL.  These  equations  are  Chen  specialized  to  the  small  signal  case, 
and  applied  to  an  experiment  we  have  recently  performed  at  Orsay. ^ 

1 .  Theoretical  development  of  the  fundamental  equations 

The  FRL  system  is  properly  described  by  the  coupled  Maxwell  and  Lorentz  force  equations, 
^^ron  these,  we  shall  derive  a  self-consistent  set  of  equations  describing  the  electron  and 
^^hc  transverse  optical  mode  dynamics.  We  use  the  dimensionless  notation  originally  devel- 
^leoed  by  W.  B.  Colson^  (In  fact  this  work  Is  s  generalization  of  Colson’s  work  to  Include 
^^ransverse  modes  and  ee  shall  stay  as  close  as  possible  to  his  original  notation). 

o 

We  begin  by  re-deriving  the  equations  for  the  slowly  varying  amplitude  and  phase  of  the 
I  ,  lOlectromagnetic  field,  but  now  including  rigorously  the  transverse  dependence.  The  field 
'"^Is  described  In  free  space  by  the  paraxial  wave  equation 
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This  equation  Is  derived  from  the  wave  equation 


0  in  the  slowly  varying 


amplitude  and  phase  approximation.  The  general  solution  of  Kq.  (1)  can  be  expressed  as  a 
linear  combination  of  a  complete  set  of  orthogonal  modes.  If  we  define  these  modes  by  the 
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complex  amplitude  K  e  where  E  Is  real  and  m  is  the  generalized  index  of  the  mode  (In  the 
two  dimensional  transverse  space"'we  consider,  m  represents  two  integer  numbers).  The  ortho¬ 
gonality  relation  reads 
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vihere  «c  have  chosen  a  convenient  nornaXlzation  which  makes  the  F.  dimensionless.  Although 
we  will  use  the  cylindrical  modes  in  the  examples,  we  proceed  witN  the  general  theoretical 
development  which  makes  no  assumptions  on  the  specific  form  of  the  modes. 

In  the  FRL,  the  coefficients  c  in  the  mode  expansion  become  tine  dependent.  Ve  wish  to 
calculate  the  evolution  of  the  amplitude  and  phase  of  these  mode  coefficients. 
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wli«T>>  we  have  made  the  definitions 
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and  where 


;(t)  -  (k  ♦  k^)z{t)  -  wt 
v(t)  -  L[(k  ♦  k^j)tf^(t)  -  k] 


dimensionless 
election  phase 

dimensionless 
resonance  parameter 

dimensionless 
interaction  time 


Mere  we  consider  an  N  period  helical  undulator  of  length  L,  magnetic  period  i  ■  2i /k  , 
peak  r.iagnetic„f Icld  B,  and  deflection  parameter  K  >  93. -I  B(gauss)i  (cm).  An  electron  beSm 
of  energy  yme^,  and  number  density  n  travels  along  the  axis  of  the”undulator ;  an  Individual 
electron  has  longitudinal  coordinate  z(t)  and  longitudinal  velocity  eg  (t)  at  tine  t.  An 
arbitrary  helically  polarized  wave  of  wavelength  i  >  2t/k,  frequency  wT  nnd  electric  field 
f  (z,t)  ■  K(z,t)e^ Interacts  with  the  electrons,  and  c  are  the  mode  coeffi¬ 
cients  in  the  expansion.  In  Kq.  (9),  <>  is  the  average  over  the'"initial  phase  c  and 

'o’'©  ° 

resonance  parameter  v^  of  the  electron  population  at  the  position  z. 

Kqs.  (3)  and  (4)  are  derived  dirsctly  from  the  Lorentz  force  equation  and  describe  the 
effect  of  the  radiation  field  on  the  electrons.  The  work  done  by  the  longitudinal  field  on 
the  electrons  is  neglected  here,  which  is  a  good  approximation  provided  that  the  modes  are 
not  too  divergent'  i/e  <<  2v^KN/y.  Kq.(5)  is  derived  from  the  Maxwell  equations  and 
describes  the  effect  or  the  electron  on  the  'radiation  field.  The  set  (3)-(9)  is  self- 
consistent.  Indeed  those  equations  are  very  close  to  being  the  most  general  classical 
equations  describing  the  PKL  dynamics.  They  apply  to  high  and  low  gain  devices  (r'  >>1  or 
r'<<l),  high  field  and  low  field  cases  (a'  >>1  or  a'  <<1),  and  Include  the  effects  of  multi¬ 
ple  longitudinal  modes  (laser  lethargy  effects)  through  the  z  dependence  of  r'  ,  E,  and, 9. 
•Slight  modifications  allow  their  extension  to  the  cases  of  a)  the  planar  undulator ^ ^ 
b)  the  tapered  undulator, c)  the  optical  klystron,'  and  d)  space  charge  effects  For 

simplicity  in  the  following  development,  we  drop  the  explicit  z  dependence  which  has  been 
thoroughly  discussed  by  Colson, and  concentrate  on  the  transverse  phenomena. 

Equation  (9>- describes  the  growth  or  decay  of  the  radiation  field  due  to  Its  interaction 
with  the  electrons.  It  shows  clearly  the  fact  that  the  growth  in  the  m*"  node  amplitude  ant 
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phase  Is  ^iven  by  the  overlap  integral  of  the  in-phase  and  out-of-phase  components  of  the 
charge  density  with  the  complex  conjugate  of  that  node,  as  one  would  expect.  We  note  that 

!<>„ 

Che  only  assumptions  made  on  the  modes  E  e  '  used  in  (3)-(b)  are  orthogonality  and  com- 
pletnness.  This  means  these  equations  are  also  valid  for  the  cases  of  waveguide  modes  and 
dielectrically  loaded  cavities. 


The  Eqs.  (3)-(5)  can  be  Integrated  numerically  to  find  the  evolution  of  the  optical  wave 
in  any  Compton  regime  FKL.  In  a  high  field  experiment,  Eqs.  (4)  and  (5)  are  nonlinear  in  a, 
and  the  wave  evolution  can  only  be  obtained  numerically.  In  this  case,  (3)-(5)  provide  a 
precise  and  efficient  technique  for  solving  the  general  problem.  In  a  low  field  situation 
such  as  we  find  at  Orsay,  however,  the  problem  becomes  linear,  and  can  be  solved  analytic¬ 
ally.  We  proceed  with  the  low  field  case  in  the  next  section. 

2.  Tho  low  fields  solution 

The  low  field  case  is  defined  by  |*.| <<  1  for  every  mode.  In  other  words,  the  electrons 
do  not  become  overbunched.  Experiments  which  operate  in  this  domain  include  the  low  field 
amplifier  experiments,  and  storage  ring  PEL  oscillators  which  saturate  by  mechanisms  other 
than  ovcrbunch ing .  The  ignition  of  any  E'EI.  oscillator  also  occurs  in  this  domain. 

Equations  (3)-(5)  can  be  solved  by  integrating  (3)  and  (4)  to  lowest  order  in  the  fields 
a  and  inserting  the  result  for  c  into  Eq.  (5).  If  the  electrons  are  uniformly  distributed 
initially  in  phase,  we  find 
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iquation  (12)  describes  a  linear  evolution  of  the  mode  amplitudes,  and  upon  integration, 
gives  the  relation 
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where  1  is  the  Identity  matrix,  and  r,.  which  is  generally  not  Hermltlan,  has  elements 
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(15) 
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The  higher  order  terms  in  g  are  proportional  to  r  and  higher  powers  of  r,  and  are  negli¬ 
gible  in  the  low  gain  case.”" 


Evidently  this  matrix  is  of  great  interest  since  multiple  passes  of  the  electron  beam 
will  result  in  multiple  products  of  this  matrix,  greatly  simplifying  the  calculation  of  tho 
modes’  growth.  We  shall  discuss  the  consequences  for  an  oscillation  experiment  in  Sect. 
2.2, 


Of  course,  one  must  keep  in  mind  that  energy  is  radiated  into  other  modes,  and  that  cross 
terns  will  mix  a  multiple  node  input.  If  the  input  beam  is  truly  monowode ,  the  power  radi- 
lated  into  Che  n^"  node  is  lower  than  that  into  the  n^  mode  by  the  ratio  |g  I ^/2Re(g  ) 
which  is  small  for  low  gain  (r  <<  1)  systems.  It  is  only  In  this  case  that  Tt  makes  sSfise 
to  speak  of  the  gain  of  a  mode.  In  high  gain  systems,  however,  the  off-diagonal  terms  can 
lead  to  substantial  enissien  of  energy  Into  the  higher  order  transverse  nodes.  If  the  Input 
bean  is  multimode,  of  course,  mode  mixing  occurs  at  all  power  levels. 


Let  us  culuulate  in  the  sinple  case  of  experlnental  interest  where  the  electron  beam 
IS  cylindrical,  and  a  good  choice  of  nodes  is  the  cylindrical  cavity  ci(;ennodes.  We  restrict 
ourselves  to  the  weakly  divergink  case  »w  ^  >>  iL  where  the  gain  takes  on  its  most  familiar 
furn.  The  node  amplitudes  and  phases  in  (13)  become  independent  of  t,  and  we  can  integrate 
the  first  tern  in  (15)  to  find  the  gain.  For  small  spread  in  resonance  parameter  over  the 
beam  pliasc  space,  the  integral  gives  the  well  known  gain  spectrum 
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In  thf  usual  exporiuental  case  (unfortunately),  lg_„|  <<  1  and  the  energy  gain  i;_  on  the 

m..de  m  becomes  ' 


This  is  exactly  the  gain  one  calculates  by  using  the  filling  factor  obtained  by  integrating 
the  ixode  profile  overlap  with  the  gain  profile,  iipeclalizing  to  the  TRM  case  with  a 
('•cussian  i.-lectron  bean  of  width  a,  we  find 


(18) 


complete  with  the  familiar  filling  factor. 

The  „  dependence  of  C  is  the  well  Known  spectral  dependence.  The  imaginary  part  of  g„ 
IS  not  nSw.  It  describes^the  phase  shift  of  the  radiation  field  as  described  by  ColsonJo 
The  ••ffect  of  the  divergence  of  the  beam  on  the  diagonal  terms  in  C  is,  to  first  order,  and 
tor  a'  filamentary  electron  bean,  the  addition  of  a  time-varying  phase  which  shifts  the 
resonance' curve  in  E<|.  (18)  by  a  constant  depending  on  the  node 
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This  means  that  the  gain  curves  of  the  modes  are  shifted  with  respect  to  each  other.  This 
phenomenon  has  been  calculated  for  the  fundamental  TFM  mode  loathe  energy  loss  approxima¬ 
tion, and  has  recently  been  observed  experimentally  at  (irsay.  It  should  be  noted  that 
for  many  practical  situations  where  the  cavity  is  optimized  for  gain  on  the  TRM  mode, 
this  expression  is  valid  for  only  the  lowest  order  mode.  The  higher  modes  become  distorted 
in  form  as  well  us  simply  shifted  in  resonance  parameter  by  equation  (19). 

In  another  paper,  we  discuss  the  optimization  of  the  cavity  in  a  low  field  oscillator.^ 

3.  Application  to  gain-vs-aperture  experiment 
3. 1  Ifescriptton  of  the  experiment 

The  gain  of  the  Orsay  FEL  has  recently  been  measured  with  the  optical  klystron  In  place^^ 
in  an  amplifier  experiment  using  an  external  argon  ion  laser  to, provide  the  coherent  mode. 

A  detailed  description  of  the  apparatus  can  be  found  elsewhere.*^  The  laser  beara  is  ana¬ 
lyzed  at  a  distance  u  from  the  optical  Klystron  after  passing  through  an  adjustable  colli¬ 
mating  iris,  which  is  centered  on  the  laser  mode  emerging  from  the  interaction  region.  The 
gain  is  measured  as  the  ratio  of  the  power  detected  in  phase  with  both  the  electron  repeti¬ 
tion  frequency  and  the  chopper  frequency  (the  amplified  power)  divided  by  the  power  in  phase 
with  the  chopper  alone  (the  incident  laser  power).  Calibration  is  performed  as  in  our  pre¬ 
vious  work.'^ 


The  tsain  is  recorded  as  a  function  of  the  iris  aperture,  and  large  variations  are  observ¬ 
ed.  One  set  of  data  points  is  reproduced  in  Fig.  1,  where  the  gain  Is  nornalized  to  its 
value  for  the  iris  coiapletely  open,  and  the  iris  diameter  is  nornalized  to  the  measured  beam 
liui.st  at  the  iris.  The  data  is  taken  at  naxinun  gain,  which  neans  v^  at  0  for  the  optical 
klystron,  and  the  laser  beam  was  carefully  aligned  to  within  about  .05  nn  of  the  axis  of  the 
oloctron  bean.  The  change  in  the  neasured  gain  as  the  iris  is  closed  means  that  the  laser 
is  not  uniformly  amplified  in  its  transverse  profile.  In  fact,  this  experiment  provides  a 
very  sensitive  technigue  for  measuring  the  power  emitted  into  the  higher  order  modes  even  in 
the  small  gain  limit  and  for  a  monomodc  input  beam.  Clearly  a  calculation  of  the  g  is 
necessary  in  order  to  explain  these  results.  In  the  next  section,  we  apply  the  theory  we 
have  developed  to  the  case  at  hand,  and  in  Sect.  3.3,  precise  comparison  is  made  between  the 
experimental  and  the  theoretical  results. 

3.3  Multinode  emission  in  a  sinclo  node  amplifier  experiment 

In  CItls  section,  we  assume  the  incident  wave  is  a  single  node  TKil  bean  with  a  weak 
tii-ld  ( |a  I  <  1),  and  perfectly  aligned  onto  the  electron  beam.  As  discussed  previously,  we 
take  the  cylindrical  eigenmodes  based  on  the  form  of  the  input  laser  bean.  Using  the  nota¬ 
tion  ot  Sect.  1,  the  input  laser  field  reaos 

i;'(T)  -  (20) 

whiT«>  the  subscript  0  refers  to  the  TKM^^  node.  From  (14),  the  output  field  K®(f)  becomes 
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.v.sMimlng  low  gain,  the  output  power  passing  through  the  ins  aperture  is 


*i'i  n  Jiif',  •Mivors  the  Iris  aperture.  The  gain  is  therefore 


/d.SK2 

J  o 


Kor  purely  cyllndrlcAl  t  •  0  modes,  G  esn  be  written 
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where  z  is  the  Rayleigh  range  qt  the  laser  mode,  z  is  the  distance  between  the  iris  and  the 
laser  beam  waist,  and  X  -  where  r^  is  the  iris  diameter.  There  are  two  interest¬ 


ing  limiting  cases: 


<:  -  2He(g^^) 


X  ♦  •  (iris  open) 


«  -  2Hc  (g„„  ♦ 


i2ptan~  — 


X  ♦  0  (iris  closed) 


It  is  obvious  from  Ki|S.  (24)-(26)  that  the  gain  changes  with  iris  diameter  in  a  way  which 
depends  on  the  magnitudes  of  the  off-diagonal  terms  in  the  gain  matrix. 

The  generalization  is  straightforward  to  the  case  or  the  multimode  input  beam,  and  to 
imperfect  alignment  of  the  laser  and  electron  beams,  although  the  calculation  becomes  more 
difficult.  This  calculation  also  applies  to  high  power  input  laser  beams  (a  >>  1)  provided 
one  keeps  in  mind  that  the  gp^  are  functions  of  a^.  ° 


GAIN 

GAIN  (iris  open) 


-- -1. 


mode  I 


Iris  diameter 
W(z) 


FiBurg  1,  The  measured  gain  as  a  function  of  the  iris  dla- 
2?*®*’  normalized  to  the  measured  beam  waist  at  the  iris. 
The  solid  points  were  taken  closing  the  iris  and  the  open 
opening  it.  The  error  bars  are  the  one  sigma 
otatistical  errors.  All  points  have  the  same  horizontal 
error  bar  which  is  shown  for  the  point  at  2.7.  The  solid 
curve  is  calculated  using  the  measured  values  for  the  elec* 
tron  and  laser  beam  sizes.  The  effect  of  each  higher  order 
node  Is  shown  by  the  dashed  curves* 


^  ^  Application  to  the  Orsay  expertipent 

ditiSns"*****^^”®"**^  ^  following  approximate 


laser 
beam : 


electron 

beam: 


-  measured  beam  waist  w  »  .67  nn 

-  wavelength  -  5145  A  ® 

-  measured  beam  waist  at  iris  w(z)  >  2.7  nm 

-  distance  from  optical  klystron  to  iris  d  •  11.6  n 

-  gaussian  and  cylindrical  with  o  »  .32  nm 


optical 

klystron:  -  Nd  ■  BO^^'*^’*** 

-  resonance  parameter  corresponding  to  naxinun  gain  with  iris  open 


Tlio  Uark  i;urve  lias  been  calculated  using  Kqs.  (12),  (13)  and  (2>1)  for  the  planar  config¬ 
uration^ taking  into  account  the  10  lowest  order  i  •  0  inodes.  The  dashed  curves  show  the 
contribution  of  each  individual  node.  The  agreenent  between  the  experinental  results  and 
the  theory  is  excellent. 

Very  sinilar  curves, were  calculated  for  other  resonance  paraneters  indicating  that  as 
expected,  the  diffraction  effects  do  not  change  much  as  a  function  of  detuning  paraneter  for 
nodes  with  low  divergence.  .  The  effects  of  the  transverse  size  of  the  electron  bean  and  the 
distance  of  the  analyzing  aperture  fron  the  FRt  are  discussed  elsewhere.  ' 
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